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Abstract 

We consider the problem of strong existence and uniqueness of a Brownian motion 
forced to stay in the quadrant by an electrostatic repulsion from the sides that works 
obliquely. When the repulsion works in a normal direction, the question has previously 
been solved by a convex analysis method. To construct the solution, we start from the 
normal case and then we use as main tool a comparison lemma. The results are remi- 
niscent of the study of a Brownian motion with oblique reflection in a wedge. Actually, 
the same skew symmetry condition is involved when looking for a stationary distribution 
in product form. The terms of the product are now gamma distributions in place of 
exponential ones. 

1 Introduction 

In the late seventies the study of heavy traffic limits in open multi-station queueing networks 
has put the question of existence and properties of the Brownian motion obliquely reflected 
on the sides of a wedge and more generally on the faces of a polyhedron. In the following 
decade there was an extensive literary output on that topic, among which we mention the 
remarkable works of Harrison, Reiman, Williams and co-authors ( [7] . [5] . |15| . |17j . [9] . [5] . to cite 
a few of them). But there is another way than normal or oblique reflection to prevent 
Brownian motion from overstepping a linear barrier. We may add as drift term the gradient 
of a concave function that explodes in the neighborhood of the faces of the polyhedron. To 
be more specific, let nj.,... , ilk be unit vectors in M. d and be real numbers. The 

state space S is defined by 

S := {x € M. d : n r .x > & r , r = 1, ...,&} . 

Let 4>i, . . . , 4>k be k convex C 1 functions on (0, oo) with 4> r (0+) = +oo for any r = 1, . . . , k. 
The potential function $ on S is defined by 



3>(x) := 0r(n r .x — b r ) . 



r=l 



Prom the results in [T] , [2] , [3] , |12] we know there exists a unique strong solution living in S to 
the equation 

dX t = dB t -V$(X t )dt 

Xq € S 1 ' 
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where B is a Brownian motion in M. d . Since 

k 

V$(x) =^4(n r .x-6 r )n rj 

r=l 

the repulsion from the faces F r = {x £ S : n r x = i r } points in a normal direction into the 
interior of S. If we introduce a set of vectors qi, . . . , q& with q^.n,. = for r = 1, . . . , k, we 
may consider a new drift function 

fc 

-5^#(n r .x-& r )(n r + q r ). (2) 

r=l 

This is a singular drift but convex analysis cannot be used as in pQ or [2] to get strong 
existence and uniqueness of the associated stochastic differential system. 

In this paper we will only take the quadrant in IR 2 as state space and an electrostatic 
drift term. We obtain strong existence and uniqueness for a large set of parameters. Our 
results are reminiscent of the thorough study in [15] about oblique reflection in a wedge. 
However a key tool in this work was an appropriate harmonic function that made a weak 
approach possible and fruitful. Thus full results were obtained, while our strong approach 
merely provides a partiel answer to the crucial question of hitting the corner. 

For theoretical as well as practical reasons, a great deal of interest was taken in the 
question of existence and computation of the invariant measure of the Brownian motion with 
a constant drift vector and oblique reflection ([8],[T7j,[9]). Under the assumption that the 
directions of reflection satisfy a skew symmetry condition, it was proved that the invariant 
measure has exponential product form density. Motivated by the so-called Atlas model of 
equity markets, some authors (|14j.[10].[llj) have recently studied Brownian motions on the 
line with rank dependent local characteristics. As this model is strongly related to Brownian 
motion reflected in polyhedral domains, the above mentioned product form of the stationary 
distribution happens to be of paramount importance. Related interesting works are also [6] 
and [4j. 

A neighboring way has been recently explored in [13]. Here the process is a Brownian 
motion with a drift term that is continuous and depends obliquely, via a regular potential 
function, on the position of the process relative to a polyhedral domain. Under the same 
skew symmetry condition as in [9], the invariant density has an explicit product form again. 
In the last section of our paper, we consider a Brownian motion with a constant drift living in 
the quadrant with oblique electrostatic repulsion from the sides. Under the skew symmetry 
condition, there is still an invariant measure in product form. Now the terms of the product 
are two gamma distributions with explicit parameters. 



2 The setting 

The general state space is the quadrant S = R+ x R + . The corner = (0,0) will play a 
crucial role and in some cases it will be necessary to restrict the state space to S° = S \ {0}. 

Let (Bt,Ct) be a Brownian motion in the plane starting from 0, adapted to a filtration 
T = [Tt] with usual conditions. Let a,(3,j,5 be four real constants with a > 0,5 > 0. We 
say that an .F-adapted continuous process (X, Y) with values in S is a Brownian motion with 
electrostatic oblique repulsion from the sides if for any t > 

X t = X + B t + af^ + (3f^ > 

Yt = n + Q + Tft + ^l ^ 
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where Xq and lo are non-negative J^-measurable random variables. Each coordinate Xt or 
Yt may vanish, so to make sense we must have a.s. for any t > 

Jo l {Xs =o}ds = f* l {Ys =o}ds = 

Jo l{x.>o}x: < 00 Jo 1 m>o}F7 < 00 • 
The drift in © is of type © with d = fc = 2 and 

4>\ (x) = -a log (x) <j) 2 (y) = -<51og(y) 
6i = 6 2 = 

ni = (1,0) n 2 = (0,l) 

qi = (o,i) qa = (f,0) 

The case where /3 = 7 = is a particular case of ([1]). In the sequel, we will note (U, V) the 
solution of the system 

U t = X + B t + af^ > 

V t = Y + Ct + 6fi% > 0. { ' 

The processes ?7 and V are independent Bessel processes (if Xq and Yb are independent 
variables). Actually, U is a Bessel process with index a — |, and the point is intanstaneously 
reflecting for C7 if a < \ and polar if a > ^. Moreover, U 2 + y 2 is the square of a Bessel 
process with index a + S, and so the corner is polar for (U, V) in any case. 

Comparison between X and U, Y and V will play a key role in the construction of the 
solution (X, Y) and the study of its behavior near the sides of the quadrant. The following 
simple lemma will be of constant use. 

Lemma 1 . For T > ; a > 0, let x\ and x 2 be nonnegative continuous solutions on [0, T] 
of the equations 

Xl (t) = Vl (t) + af^ 
x 2 (t) = vtW+aJt^fo 

where v\, v 2 are continuous functions such that < Vi(0) < v 2 (0), and v 2 — v\ is nondecreas- 
ing. Then x\{t) < x 2 (t) on [0, T\. 

Proof. Assume there exists t G (0, T] such that x 2 (t) < x\(t). Set 

r := max{s < t : x\(s) < x 2 (s)} . 

Then, 

x 2 (t) - xi(t) = x 2 (t) - xi (r) + (v 2 (t) - vi(t)) - Mr) - t*(r)) + / T '(^ - ^)ds 

> 0, 

a contradiction. ■ 
We will also need the following form of the results in [lj,[2j,[3j,|12j on multivalued stochas- 
tic differential systems. 

Proposition 2 . Let a > 0, 5 > 0, a = (o~y,i,j = 1,2) a 2 x 2-matrix, (Bi,B 2 ) a Brownian 
motion in the plane, b\ and b 2 two Lipschitz functions on R 2 , Zq and two J-q -measurable 
nonnegative random variables. There exists a unique solution (Z 1 ,^ 2 ) to the system 

Z\ = zl + olBl+elBt + a^jk + tib^ZlZ^ds 

Z 2 = Zl + alBj+alB^ + S^^ + ^ZlZ^ds W 

with the conditions Z\ > if 5 = and Z\ > 0, Z 2 > if 5 > 0. 
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3 Avoiding the corner 

We shall see in the next section that existence and uniqueness to the solution to ([3]) are easily 
obtained as soon as the solution process keeps away from the corner. Thus the question of 
hitting the corner in finite time is of great interest. 

Theorem 3 . Let (X, Y) be a solution to |3p in the interval [0, r] n [0, oo) where r is a 
T -stopping time. We set 

t° : = inf{t G (0, t] n (0, oo) : (X t ,Y t ) = 0} . 

Then P(r° < oo) = if one of the following conditions is satisfied: 

1. d:fi>0andj>0 

2. C 2a :a>\ and /3 > 

3. C 2b : S > \ and 7 > 

4- C3 : There exist A > and fx > such that 

• Xa + /i7 > 

• A/3 + fi5 > 

• A(Aa + ^7) + n(X/3 + nS) - ±(A 2 + fi 2 ) > -2 v / A/u(A/3 + fj,5)(Xa + ^7). 

Proof. 

Condition C\. From Lemma 1 we get > Ut, Yt > Vj. where (U, V) is the solution to 
and we know that is polar for (U, V). 

Condition C 2a (resp. C 2 b)- Prom Lemma 1 we get Xt > Ut (resp. Yt > Vt) and in this case 
is polar for U (resp. V), so U t > (resp. V t > 0) for t G (0, r] n (0, 00). 
Condition C3. For e > let 

a e = l {( x ,y )=o} inf{* G (0, r] n (0, 00) : X t + F t > e} 
r°< e = inf{t G (<T e , r] n (0, 00) : {X t , Y t ) = 0} . 

As e I 0, cj e I on {r > 0} and r°' e It . We set S t = XX t + ^ for t G [0, r] D [0, 00). From 
Ito formula we get for t G [a e , r 0,e ) n [0, r] 

log 5* 

= logS^ + £ Mto^ + (Aa + M7) £ ^ + (A/3 + ^) £ ^ - I(A 2 + /, 2 ) & Jf 
= logS^+M. + ^^fd, 

where M is a continuous local martingale and P(x, y) is the second degree homogeneous 
polynomial 

P(x, y) = A(A/3 + ii5)x 2 + /x(Aa + ^)y 2 + (A(Aa + ^7) + p(\p + ^J) - ^(A 2 + ^ 2 ))xy . 
Condition C3 is exactly the condition for P being nonnegative in S. Therefore 
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and so 



, ' e/Xr P(X S ,Y S ) 



As t — > T°' e A r, the local martingale M either converges to a finite limit or oscillates between 
+00 and —00. It cannot converge to —00 and thus S T o,e AT > on {r°' e A r < 00}, proving 
that P(r°' e < 00) = and finally P(t° < 00) = 0. ■ 
Example. When a = 5 and \0\ = I7I, Condition C3 is satisfied (with A = /i) if 

• 2 < a - \ when = -7 , , 

• -0<a-\ when /? = 7 < . U 

We may also be interested in hitting a single side. Then we set 
r° := inf {t G (0, r] n (0, 00) : X t = 0} 

r£ := inf {t G (0, r] n (0, 00) : Y t = 0} . l j 

We already know that P(t^ < 00) = if a > \ and /3 > 0. In a symmetric way we can prove 
that P(t^ < 00) > if a < \ and < (= 1 if moreover r = 00 a.s.). If we know that the 
corner is not hit and a > |, we can get rid of the nonnegativity assumption on 0. 

Proposition 4 . Assume P(r° < 00) = 0. Ifa>\, then P(r£ < 00) = 0. 

Proof. For rj > let 

4 = l{x =0} inf {t E (0, r] n (0, 00) : X t > rj} 
= inf{t G {6 V X , t] n (0, 00) : X t = 0} . 

As ?? I 0, 0£ I on {r > 0} and t% v | t° . For i G [0^, tJ") n [0, r], 

logX, = logX e ,+4^ + (a-i)/4^+^4 A. (8) 

Since P(r° < 00) = 0, we have Y T o, v > on {t^ 1 < 00}. On this set, 

0,1) 0,1) 

r * ds Z" 1 "^ (is 

— < OO , / — < OO 

3^7 0,V 



and X s > on [0jp r x ?? )> wn i cri proves that 

a r x ds 

- T1 ->-oo. 

As t — > Ty 7? , the local martingale in the R.H.S. of ([8]) cannot converge to —00. This entails 
that P(r^ ?? < 00) = and therefore P(r£ < 00) = 0. ■ 

4 Existence and uniqueness 

We now proceed to the question of existence and uniqueness of a global solution to ([3]). We 
consider separately the three cases: > and 7 > 0, > and 7 < 0, < and 7 < 0. 



5 



4.1 Case /3 > and 7 > 
This is exactly Condition C\. 

Theorem 5 . Assume (3 > and 7 > 0. 

1. There is a unique solution to in S°. 

2. There is a solution to §3$ in S starting from 0. 

3. If a5 > /3j, there is a unique solution to |5}) in S. 
Proof. 1. Let a > 0, e > and define for (x, z) G W + x R 

ifj e (x, z) : = max( — - — , — ) . 

73; + z ae 

This is a Lipschitz function. From Proposition [2] we know that the system 

Z\ = -7X0 + a(Y + l{y 0=0 }a) - jB t + aC t + a(a5 - ^7) J V>eP^, Z e s )ds 
has a unique solution. Let 

rf := inf{* > : jX% + Z\ < ae} . 

If < 7] < e < a we deduce from the uniqueness that (X e , Z e ) and (X 11 , Z 11 ) are identical on 
[0, Ty ] . Patching together we can set 

X t := lim e ^o -X* 

Y t := lim e _ I( 7 ^ + Z|) 



on {Y > 0} x [0,r y ), where 



t y := limrf . 

e-*-0 



On this set, (X,Y) is the unique solution to ((3|). As we noted in the proof of Condition Ci, 
we have X t > U t and If > V t . Therefore, on {Y > 0} n {r y < 00}, 

tf Y %<tfv;<°° and Jo 4 g</o r "g<- 

and we can define 



X T o := lim^oX, = X + S r o + a + /3 § 

y T o := lim^y* = ^o + C r o +7/o T "* + ^o 4 ft- 



(10) 



We have Y^o = and as is polar for (U, V), then X T o > 0. In exactly the same way we can 
construct a solution on {Yq > 0} in the interval [Ti,T2], where T\ = t y , I2 = inf{i > T% : 
X t = 0}. Iterating, we get a solution on {Yq > 0} n [0, hin^^oo T„) where 

T 2p := inf {* > T 2p _! : X t = 0} 
T 2p+1 := inf{t > T 2p : Y t = 0} . 

On {Yq > 0} n {lim^ooTn < 00} we set X limrwooTn := ]hn v ^ 00 X T2p = and F^^^ := 
limp^oo Yr 2p+1 = 0. The polarity of entails this is not possible in finite time and thus 
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linin-j.oo T n = oo. So we have obtained a unique global solution on {Yq > 0}. In the same 
way we obtain a unique global solution on {Xq > 0} and as P((Xo, Yq) = 0) = the proof is 
complete. 

2. Assume now Xq = Yq = 0. Let (y n )n>i be a sequence of real numbers (strictly) 
decreasing to 0. Consider for any n > 1 the unique solution (X n , Y n ) to the system 

x t = Bt + a Jo Jfr + /? J * ^ 

From the continuity of It there exists a positive stopping time p such that a.s. l^ n+1 < Y™ 
on the interval [0,p]. From Lemma [J we deduce that X™ +1 > X™ on [0,p]. Let 

r := inf {t > : X™ +1 < X™} . 

Using again Lemma CD we obtain Y t n+1 < Y t n on [0,r]. On {Y r n+1 = Y r n } n {t < oo}, since 
X™ +1 = X™, we note that (X™, Y T n ) € S" . Therefore the uniqueness of the solution starting 
at time r implies that X™ +1 = X t n and Y™ +1 = F t n on [r, oo). On {Y r n+1 < Y™} n {r < oo}, 
the same proof as above entails there exists p' > such that X™ +1 > X™ on [r, r + p'], a 
contradiction to the definition of r. Therefore X™ +1 > X™ and < Y t n for any t G [0, oo). 

We define 

X t := lim t X? Y t := hm 1 Y t n . 

n— >oo n— >oo 

As Y t n > V t where (U,V) is the solution to (@J) with Xq = F = 0, we have 
Xt = -Bt + a lim^oo f Q ^c+P linin->oo J 4^ 

< oo 

and also 

Y t = lim^oo y n + C t + 7 lim^oo / Q ^ + 5 lim^oo J Q ^ 
< oo . 

3. Assume finally a<5 — /?7 > 0. As the conclusion holds true if (3 = 7 = 0, we may 
also assume that /3 > 0. Let (X, Y) be the solution to ([3]) with Xo = Yq = obtained in 
the previous paragraph and let (X',Y') be another solution. Considering (X n ,Y n ) again 
and replacing (X n+1 ,Y n+1 ) with (X',Y'), the previous proof works and we finally obtain 
X' t > X t and Y{ < Y t . Then, 

(S(X t - X't) - (5(Yt - Y()f 
= 2fl{5(X s -X> s )-P{Y s -Y>)){ a 5-M){jc s -jtr)ds (11) 
< 

and thus X[ = X t and Y( = Y t , proving uniqueness. Same conclusion with a small change in 
equation (fTT]) if 7 > 0. ■ 

4.2 Case > and 7 < 

Theorem 6 . Assume (3 > 0, 7 < and one of the Conditions C2 a or C3 is satisfied. Then, 
there exists a unique solution to |3j) in S°. 
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Proof. The proof is similar to the proof of 1 in Theorem [5j The only change is that now 
Y t < V t . Therefore, on {Y > 0} n {r£ < oo}, 



_<^_r -c r „- 7 y o -<oc 

and we can define and Y T o as previously. ■ 
4.3 Case < and 7 < 

In this case we can give a full answer to the question of existence and uniqueness. Our 
condition of existence is exactly the condition found in [16] for the reflected Brownian in a 
wedge being a semimartingale, i.e. there is a convex combination of the directions of reflection 
that points into the wedge. 

Theorem 7 . Assume (3 < and 7 < 0. 

1. If aS > f3j, there exists a unique solution to ([3p in S. 

2. If a5 < /?7, there does not exist any solution. 

Proof. 1. Existence. Let (h n ,n > 1) be a (strictly) increasing sequence of bounded positive 
nonincreasing Lipschitz functions converging to 1/x on (0, 00) and to +00 on (— 00, 0]. For 
instance we can take 

h n (x) = (l-^)i on [loo) 

= w-l on (-00, ±] . 

We consider for each n > 1 the system 

X" = X + B t + af*^ + Pf*h n (Y?)ds 

Y? = Yo + Ct + jfthnlx^ds + Sj*^. 1 ' 

From Proposition [2] we deduce there exists a unique solution to this system. As h n+ i(Xo) > 
h n (Xo) and h n +i(Yo) > h n (Yo), by continuity there exists an a.s. positive stopping time 
p such that h n+1 (X™ +1 ) > h n {X™) and h n+1 {Y s n+1 ) > h n (Y s n ) on [0,p]. A first double 
application of Lemma CD shows that < X t " and Y t n+1 < Y t n on the interval [0,p]. Now 

we set 

r := inf{s > : X r s l+1 > X™} . 

We have r > p > and a second application of Lemma [1] shows that Y™ +1 < ^j™ on [0, r]. 
As in the proof of 2 in Theorem [5] and since /i n +i(Y r n+1 ) > h n (Y™) on {r < 00}, a third 
application of Lemma [1] shows that there exists p' > such that X™ +1 < X" on [r, r + //], a 
contradiction to the definition of r. Thus r = 00 proving that on the whole [0, 00) we have 
X™ +1 < X t n and Y t n+1 < Y t n . We can set for any t € [0, 00) 

X t := lim X t n and Y t := lim Y t n . 

n—too n— >oo 

If a5 > /3j, there is a convex combination of the directions of repulsion that points into the 
quadrant, i.e. there exist A > and p > such that Xa + /try > and /u<5 + A/3 > 0. For 
n > 1 and £ > 0, 

AC/t + /^ > AXf + /xF t " 

> AX + pY + XBt + + (Aa + /n) £ ^ + W + ^ /o & ■ 
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Letting n — > oo in (|13p we obtain 

* ds , /•* ds 

— < oo and / — < oo . 

-^-s JO *s 

Then we may let n go to oo in f 1 1 2 [) proving that (X, Y) is a solution to ([3]). 

Uniqueness. Let (X',Y"') be another solution to ©. Replacing (X n+1 ,Y n+1 ) with (X',Y') 

we follow the above proof to obtain for t € [0, oo) and n > 1 

X;<^" and Y/<Y™ 

Letting n->oowe conclude 

X[ < X t and Y/ < Y t . 
With the same A > and fi > as above, 

(X(X t - X' t ) + n(Y t - Y{)f 
= 2 f*(\(X s - X' s ) + [i(Y s - YD) [(Aq + ^ 7 )(£ - + + A/3)(£ - £)]<fe 
< 

and therefore X^ = X t , Y( = Y t . 

2. If a5 < /?7 there exist A > and > such that Aa + /X7 < and fiS + A/3 < 0. Thus, if 
(X, Y) is a solution to (J3J), 

< AX + ^Y* < AX + fiY + \B t + \iC t . 

This is not possible since the paths of a Brownian motion are not bounded below. So there 
is no global solution. ■ 

5 Product-form stationary distribution 

We introduce an additional constant drift (— fi, —v) in the quadrant and consider the system 

with the conditions Xt > 0, Yj > 0. If now 

V = Y + Sj^-ut U5) 



with Ut > 0, Vt > 0, we can check that is still polar for (U, V) (as well, is polar for U if 
a > I and for V if <5 > i). Therefore the results of the previous sections are still valid for 
the solution to (|14p . We are now looking for conditions on the set of parameters in order to 
obtain a stationary distribution for the Markov process (X, Y) in the form of a product of 
two gamma distributions. 

Theorem 8 . Assume there exists a unique solution to |7^p in S° or in S. This process has 
an invariant distribution in the form T(a,c) (g> T(b,d) if and only if 

• a/3 + 7<5 = (skew symmetry) 

• a = 2a + 1, b = 25 + 1 

. C = 25^±P, d = 2a^f 

ad—p^f ' ad— pry 

• \xa + 1^7 > 0, + v5 > . 
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Proof. Let 

p{x, y) = x a - x e- cx y h - x e- dy for x > 0, y > . 
The infinitesimal generator of the diffusion (|14p is given by 



_ 1 , d 2 d 2 a f3 d 7 5 d 
2 dx 2 dy 2 x y dx x y dy ' 

By a density argument it is enough to check that 

00 roo 

Lf{x,y)p(x,y)dxdy = 



J 

for any f(x,y) = g{x)h{y) with g,h € C,?((0,oo)). Integrating by parts, we get 

OO /'OO /'OO /"OO 

/ L(gh)(x,y)p(x,y)dxdy = / g(x)h(y)J(x,y)dxdy 
io jo io 

where 

J( x , y ) = p( x , y) [A + Bx' 1 + Cx~ 2 + Dy" 1 + Ey~ 2 + 

with 

A = \(? + i(i 2 - pc-vd 

B = — (a — l)c + /i(a — 1) + ac + 'yd 

C = i(o- l)(a-2) -a(a-2) 

D = -(& - l)d + zy(6 - 1) + /3c + <5d 

£ = i(6_i)(6_2) -5(6-2) 

F = /3(a-l)+ 7(6-1). 

Letting J 4 = i? = C = L' = £ , = i ? = 0we obtain the skew symmetry condition a/3 + 7<5 = 
and the specified values for a, b, c, d. The last condition in the statement of the theorem is 
written out so that the invariant density p is integrable on S. Note that when p = v = 0, the 
invariant density p(x, y) = x 2a y 2S , which is not integrable, does not depend on the obliqueness 
parameter f3. ■ 

Acknowledgement. The author is grateful to N. Demni for drawing his attention to reference [13] . 
which motivated this reseach. 
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